BBGKY EQUATIONS, SELF-DIFFUSION AND 1/F NOISE 
IN A SLIGHTLY NONIDEAL GAS 



YURIY E. KUZOVLEV 

Abstract. The hypothesis of "molecular chaos" is shown to fail when applied to 
spatially inhomogeneous evolution of a low-density gas, because this hypothesis is 
incompatible with reduction of interactions of gas particles to "collisions" . The fail- 
ure of molecular chaos means existence of statistical correlations between colliding 
and closely spaced particles in configuration space. If this fact is taken into account, 
then in the collisional approximation (in the kinetic stage of gas evolution) in the 
limit of infinitely small gas parameter the Bogolyubov-Born-Green-Kirkwood-Yvon 
(BBGKY) hierarchy of equations yields an autonomous system of kinetic equations 
for the many-particle distribution functions of closely spaced particles. This system 
of equations can produce the Boltzmann equation only in the homogeneous case. It 
is used to analyze statistical properties of Brownian motion of a test gas particle. 
The analysis shows that there exist fluctuations with a 1/f spectrum in the diffusivity 
and mobility of any particle. The physical cause of these fluctuations is random- 
ness of distribution of particles' encounters over the impact parameter values and, 
consequently, randomness of the rate and efficiency of collisions. 

In essence, this is reprint of the like author's paper published in Russian in [ Zh. 
Eksp. Teor. Fiz. 94 (12), 140-156 (Dec. 1988)] and translated into English in [ Sov. 
Phys. JETP 67 (12), 2469-2477 (Dec. 1988)] twenty years ago but seemingly still 
unknown to those to whom it might be very useful. The footnotes contain presently 
added comments. 



I. Introduction 

The today's kinetic theory of weakly nonideal gases as before rests on the antiquated 
hypothesis of "molecular chaos" which asserts that the particles entering a collision are 
statistically independent and which makes it possible to reduce the exact Bogolyubov- 
Born-Green-Kirkwood-Yvon (BBGKY) equations to the classical model Boltzmann 
equation [HE], [3]. Meanwhile, the molecular chaos never has been proved and, in fact, 
can be justified only for the special case of spatially homogeneous gas evolution [21 [31 
HI El E]. As for the general case, any reasonings in favor of the molecular chaos involve 
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other arbitrary assumptions. For example, it is sometimes identified with those fact 
that in a sufficiently low-density gas the colliding particles do not have an intersecting 
dynamic past. However, as was emphasized many years ago [3], generally the absence 
of dynamic correlations by no means implies the absence of statistical correlations 
from the standpoint of the probability laws which pertain to an ensemble of dynamic 
systems. It is also incorrect to identify molecular chaos with the decoupling of statistical 
correlations for infinitely far spaced particles, since in actuality the particles arrive at 
a collision not from infinity but from a distance which is only on the order or less than 
the mean free path Ao (moreover, in practice, molecular chaos is being postulated even 
for distances on the order of the interaction radius ro <C Ao )l • 

On the other hand, it is not difficult to indicate why the molecular chaos can fail in 
inhomogeneous non-equilibrium situations. Notice that, first, in spatially inhomoge- 
neous gas configurational (spatial) dependencies of distribution functions (DF) carry 
statistical information not only about the instantaneous coordinates of gas particles but 
also, indirectly, about their past diffusive displacements, or "Brownian paths" (since 
the non-homogeneity constitutes a natural reference scale for the displacements). Sec- 
ond, the displacement of each particle is closely correlated with fluctuations in the rate 
of collisions of this particle and consequently is correlated, to the extent of duration 
of these fluctuations, with next collisions. Therefore the pair (two-particle) DF for 
pre-colliding particles (i.e. particles going into mutual collision) actually represents a 
conditional probability distribution under the conditions that a new collision realizes 

1 Another example of the arbitrary assumptions is replacement of the BBGKY hierarchy by so 
called "hard sphere BBGKY hierarchy" which has no substantiation but was exploited in the Lan- 
ford's attempt to substantiate the Boltzmann equation (although for absurdly small evolution time 
only) [ O.E.Lanford, "Time evolution of large classical systems", in "Dynamical systems, theory and 
applications", ed. J.Moser, Lectures Notes in Physics, vol.38, 1975; "On a derivation of Boltzmann 
equation", in "Nonlinear phenomena. 1. The Boltzmann equation" , eds. J.L.Lebowitz, E.W.Montroll, 
N.-H., Amsterdam, 1983; H.Spohn, "Theory of fluctuations and Boltzmann equation" , ibid.]. In fact, 
the "hard sphere BBGKY hierarchy" is a "hand-made" probabilistic model which does not follow 
from the Liouville equation even in the limit of infinitely hard repulsion. The present paper helps to 
comprehend why it is so (for detail, see [Yu. Kuzovlev, "On Brownian motion in ideal gas and related 
principles", arXiv: 0806.4157] and footnotes below). 
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and it takes place in a given small space-time region. Even because of this circumstance 
alone, generally the mentioned DF can not be factored into the product of one-particle 
DF which would furnish only unconditional information about the coordinates and 
displacements. 

The question thus is the extent to which the fluctuations in the collision rate of 
an arbitrary gas particle are "long-living". A natural source of these fluctuations is 
the randomness of the impact parameter in an encounter of two particles (we will 
be using the word "encounter" to mean both a factual interaction of particles and 
their mere passing at a relative distance comparable with the interaction radius r ). 
It is quite obvious that this source does not reduce completely to the ordinary local 
gas density fluctuations. Being dependent on a random distribution of the particles' 
encounters over values of the impact parameter, the actual rate of collisions of any 
given particle can randomly vary from one time interval to another. But, at the same 
time, a thermodynamic state of the gas is indifferent to these variations, since anyway 
they do not interfere with a local thermalization of the gas during time on order of 
the mean free path time r . Consequently, there are no relaxation mechanisms which 
would tend to establish some certain distribution (histogram) of the impact parameter 
values, and in this sense the relaxation time (or smoothing time) of this distribution 
is infinitely long. Such the reasonings show that fluctuations in the collision rate (and 
thus in the gas kinetic coefficients associated with this rate) are of a long-living "flicker" 
nature EH 13 M EE] • 

To deal with these fluctuations we must abandon a priori molecular chaos, i.e. treat 
the pair DF, for particles which are encountering one another (in the sense explained 
above), as an autonomous statistical characteristics of gas evolution, which represents 
a local ensemble-average rate of encounters and thus rate of collisions proper. Accord- 
ing to the BBGKY equations, the evolution of this DF is coupled with evolution of 
the higher-order DF for "clusters" of n > 2 relatively close encountering particles. 
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One might suggest that, taken together, they represent temporal statistics of impact 
parameters of particles'encounters and thus statistics of the rate of collisions. 

In Sec. 2 we will show that in the framework of the coarsened "collisional" (kinetic) 
description of the particles' interaction the BBGKY hierarchy generates a separate 
system of evolution equations for just mentioned special DF of particles under mutual 
encounters. At that, the structure of these equations is such that in spatially inhomo- 
geneous case it forbids Boltzmann's molecular chaos. The only possibility in general 
is a weakened version of the molecular chaos hypothesis (examined in Sec. 3) which 
incorporates the inter-particle statistical correlations in configuration space. 

Such a weakened hypothesis, however, is sufficient for deriving a closed (although 
infinite) system of kinetic equations. As is shown in Sec. 4 by the example of self- 
diffusion, these equations predict flicker (1/f) fluctuations in the transport coefficients 
of a gas. This result, discussed in Sec. 5, supports the fundamental conception of 1/f 
noise which was first proposed in [3 [8]. 

2. Collisional approximation 

Since the BBGKY hierarchy can not be solved rigorously, we naturally appeal to 
the idea, which was suggested by Bogolyubov in [I], about possibility of asymptotic 
separation of "collisional" and "kinetic" space-time scales in the low-density limit 
(Ao = const , p = r /A ~ prjj — > 0, where p is mean gas density). In other words, 
following [1], let us assume that at a sufficiently late stage of the gas evolution all 
many-particle DF F n possess, along with F\ , only slow time dependence character- 
ized by "kinetic" time scales > t = A /t>o (v = y/T/m is thermal velocity). In order 
to implement this idea into practice we have to specify the approximate asymptotic 
form in which we are seeking a solution of the BBGKY equations. For this purpose, 
Bogolyubov in pQ introduced the assumption that all of the DF F n are local, in re- 
spect to time, functionals of the one-particle DF Fx . That assumption makes a use 
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of the molecular chaos hypothesis unavoidable, although the BBGKY equations by 
themselves in no way impose this hypothesis. The "slowness" of F n , however, gives a 
longer list of possibilities. We will accordingly discuss a less rigid formulation of the 
separation of scales. 

To consider DF for closely spaced particles, it is natural to express the inter-particle 
distances q^ = q^ — qj ( ^ are coordinates) in units of r , while the position of the 
center of mass of a cluster as a whole, q^ = (1/n) Y^j=i Qj > i n units of A . Putting 
the particle velocities in dimensionless form by dividing by thermal velocity scale v q , 
and putting the time in dimensionless form by the mean free time r (in accordance 
with the presumed slowness of changes of DF), we can write the following expression 
for the volume- normalized DF: 

F n (t , 8 W, qij , VJ ) = F n (- , ^ , *i , V -l) . (2.1) 

The separation of scales presupposes that in a certain asymptotic sense the "reduced" 
DF of close-lying particles, F n , do not depend on the gas density rho , i.e. do not 
contain the gas parameter \i = pr^ as a special independent argument. Let us examine 
the consequences of this - still preliminary - suggestion by substituting (12 .ip into the 
BBGKY equations. 

To exhibit the scale r explicitly, it is convenient to specify the inter-particle in- 
teraction force to be (T/r ) f(qij/r ) . At that, we can always choose r and A in 
such a way that the relation A = holds. Let us introduce the designations 

Zn — q /Ao and Xij = qij/ro, while retaining the old notations t and Vj for new 
dimensionless time £/tq and velocities vj/vq. The BBGKY equations can then be 
easy put in the following form: 

(d d 1 \ ~ n d f ~ 

di + Un dz~ + Jt L ' n ) Fn = ^ ~dv- J f( Xn+lj ^ Fn+1 dXn + l i dVn + l = Jn ' ( 2 ' 2 ) 

where operator L' n (which acts on the functional dependence on X{j and Vj only) 
is the Liouville operator of the relative motion and interaction of n particles, which 
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results from the complete n-particle Liouville operator by eliminating the center-of- 
mass motion, and 




is the center-of-mass velocity. 

It is seen from equations (12.21) that formally strict independence of F n on // would 
imply a supplement to the equations in the form of the requirement 

L' n F n = 0, (2.3) 

which just excludes contributions of "fast" relative interactive motion of particles. The 
physical meaning of this requirement is easy understandable. It asserts that different 
dynamical states, which realize in the course of an encounter of n particles on the same 
phase trajectory in n-particle phase space, have the same probability in the statistical 
ensemble under consideration . In other words, different dynamic stages of the same 
collision (in particular, in- and out- states) are represented by the statistical ensemble 
with equal weights. 

In essence, this statistical property is obligatory attribute of such a gas evolution 
which allows coarsened description in terms of collisions (when details of geometry and 
time history of the interaction process are replaced by indication of only input and 
output of a momentary "collision"). Therefore we can expect that exact solution of 
the BBGKY hierarchy asymptotically satisfies the requirement (12. 3p . and thus (12.31) 
should be thought of as more adequate basis of the collisional approximation than that 
proposed in [Tj (indeed, the alternative to (12.31) would be unfitness of the concept of 
collision at all, which would contradict the elementary physical logics). 

At the same time, undoubtedly, the equality (12. 3j) never holds in literal rigorous 
sense, at least because the presupposed scale separation concerns only not too long 
inter-particle distances, at any case |gy| <C Ao ( \xij\ <C /U _1 ). We will thus move on 
to a more correct treatment of the separation. For this purpose it is quite sufficient 
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to understand (I2.3P as the condition that the quantity /i -1 L' n F n (or L' n F n in the 
original dimensional form) is small in comparison with the other terms of the n-th 
BBGKY equation. Furthermore, it is sufficient if (I2.3P holds only on the average over 
some region in the space qij with a linear size a much larger than tq but much 
smaller than Ao . A natural (and unambiguous in order of magnitude) choice for a if 
the average distance between neighboring particles: a = p -1 / 3 (in the dimensionless 
form, a/ro = /i~ 1//3 ). 

In the limit /i — > this region (the "collision volume") becomes infinitely large on 
the scale of r but it shrinks to a point at the scale of Ao • Then one can neglect 
the vanishingly small (< a/Ao ~ /i 2//3 ) difference between the centers of mass of the 
configurations on the left and right sides of (12.21) . and replace the chain of variables 
z n by the single common variable z : the coordinate of a physically small collision 
volume. The belonging to the same such volume will be taken below as the criterium 
of closeness of particles. 

Let us denote the mentioned averaging operation by the overline, and the result of 
the averaging of F n by A n = F n . By virtue of this definition of DF A n , any of 
A n = A n (t, z, vi, v n ) depend only on t , z and the velocities Vj and characterizes 
a local mean (ensemble-average) density of the number of n-particle encounters. At 
that, according to the aforesaid, now instead of (12. 2p we have 



(or L' n F n in the dimensional form). Due to this equality the equations (12.21) turn into 



n 







(2.4) 




(2.5) 



where V = d/dz , and (as in (12. 4p ) the limit fi — > is taken in mind. 
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It is thus clear that in general inhomogeneous case, when VA n ^ , a solution of 
equations (12.51) can not be written as the product of one-particle DF: 

n 

A n (t,z,vi, ...,v n ) ^ Yl Ai{t, z > v i) > 

since the inertial terms u n VA n constantly generate statistical correlations between 
close-lying particles, due to their joint drift (together with their collision volume) rela- 
tive to the inhomogeneity. As the consequence, the circumstance that particles belong 
to the same encounter or collision event already establishes statistical correlations be- 
tween them. At that, of course, from the probabilistic point of view, there is no 
principal difference between encounters and collisions proper. 

We thus arrive at the conclusion that in non-homogeneous situations the collisional 
approximation - by virtue of its very nature - contradicts the hypothesis of molecular 
chaos, since in the language of collisions the relative motion of colliding particles be- 
comes an inner constituent part of the collision as a whole and therefore automatically 
excluded from the equations for DF what characterize number density of collisions. 
As a result, even in the hydrodynamic stage, the evolution of gas is described by the 
infinite system of equations (12. 5p . 

But, on the other hand, the collisional character of gas evolution does not prevent 
factorization of the n-particle DF F n (t, q 1 , q n , Vi, v n ) if the particles are suffi- 
ciently far apart from each other, i.e. are not close in the above-defined sense. Indeed, 
if in the limit tq/Xq — > ( Ao = const) one keeps the inter-particle distances qij fixed 
in units of Ao , instead of r , then the BBGKY equations reduce to such the equations 
for F n which have the factored solution F n = Y[j Fi(t, Qj ,Vj) . However, in the (right 
side of) equation for F\ we see the pair DF taken at quite different type of limit, when 
qi2 is fixed in units of r , i.e. F 2 | g2 = gi , which leads to the hierarchy (I2.5p p . 



2 Notice that our reasonings nowhere appeal to details of the interaction potential (presuming only 
that it is short-range enough). Therefore nothing prevents us to extend our conclusions to the limit 
of the hard-sphere interaction. Moreover, there is no alternative, since only ansatz like (|2.3|) ensures 
that probabilities are conserved during collisions. 



The switch to the common spatial variable z in (12. 5p . of course, presumes that 
| g (n)_g(n+i)| ^ i anc j nfl 3 ^ ^3 ^ w h ere / j s characteristic scale of the non-homogeneity 
and na 3 is characteristic volume of a cluster of n close particles. In the limit under 
consideration, both these requirements are satisfied by an infinite margin if I > Ao 
(since then na 3 /I 3 ~ n/i 2 — > ) and thus do not restrict the number of DF A n which 
are "tied" to a given coordinate z . 

3. Weakened molecular chaos 

The requirements (12.41) are main tools of construction of a collisional approximation: 
due to them the right sides of (12.21) and (12. 5p can be reduced to the collisional form. 
In particular, with n = 2 , Eq J2.4l becomes 

fi- l L' n F n = CTV 1 f L' n F n dq 21 = j L' n F n dx 2l =0 

or, after we take the limit \i — > , 



/ 

J X 



(- 


--)} 


\dv 2 


dvxj _ 



F n dx 2 i = . 



|z2i|<oo 

With the help of this equality the right side J\ = J\ in the first of equations ( 12. 2ft and 
( 12. 5ft transforms to integral 



■h = J dv 2 (t> 2 - Vi) j) dsF 2 



over an infinitely remote surface |x2i| = oo (with ds being its normal vector), so that 
Ji is determined by the particle flow into the "collision volume" \x 2 i\ < oo from the 
surrounding gas. 

Depending on the sign of the scalar product (v 2 —Vi)-ds , the DF F 2 represents either 
in- or out -state of particle 2 with respect to particle 1. Let us denote by A % 2 the 
values of lim^o-^2 on that part of the boundary surface |g2i| — o, which corresponds 
to m -states. The boundary values for the out -states can then be expressed in terms 
of the A l 2 n with the help of the two-particle scattering matrix. After that, J\ acquires 
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the standard form of the collision integral: 



8 „\ 



1 — 



/ 



dv 2 S 12 A 2 r 



Here and below Sij is the ordinary "Boltzmann collision operator" for the collision of 
particles i and j . The action of this operator is defined by [U [2j [3] 



velocities which correspond to the final velocities Vi and Vj . 

Analogously, we can use (12.41) with n > 2 to perform similar transformations of the 
integrals J n . In the limit jj, — > the functions A n are determined by the average 
of F n over an infinite (3(n — l)-dimensional) region of dimensionless inter-particle 
distances. Therefore the result of this averaging represents only such (pre- or post- 
collisional) configurations of n particles where none of them are just now in a collision. 
Correspondingly, only two-particle collisions (between some of n lefthanded particles in 
(12. 5p and "external" righthanded (n + l)-st particle from the rest of the gas) contribute 
to J n . We thus find what could have predicted earlier: 



where A™ +1 is the boundary DF (similar to A 2 n ) representing configurations with the 
external particle which always is in infinitely remote in -state in respect to other n 
particles. To underline the particular role of the external particle, we will distinguish its 
velocity among other arguments of A™ +1 and write A™ +1 = A™ +1 (t, z, v\ , v n \ v n+ \) . 

In order to transform ( 13.1 p into a closed system of equations, we have to relate the 
right-side boundary DF A l ™ +1 to the left-side functions. In this step - after transition 
to the collision integrals - we need to invoke the concept of molecular chaos. Concretely, 
let us assume that the external (n+ l)-st particle, due to its just noted specificity, has 




where b is the two-dimensional impact parameter vector, and v • and v'j are the initial 




(3.1) 
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no velocity correlations with the other particles: 



A™ +1 (t, Z,V! , ...,v n \ v n+ i) = Ai(t, z,v n+1 ) A' n (t, Z,Vx , ...,v n ) . 



However, this velocity factorization does not mean absolute statistical independence, 
since it still allows a spatial correlation, by virtue of which the function A' n may 
differ from A n (according to the definition, A' n is the conditional n-particle DF 
corresponding to the condition that a collision with an additional particle takes place). 

In the "pure" form the correlation of particles in the configuration space is described 
by the DF integrated over all velocities. Since in all the configurations under consid- 
eration the particles are infinitely close together from the standpoint of the scale Ao , 
the degree of their spatial correlations in all these configurations should be the same. 
This statement is expressed by the equality 



In essence, it claims the conservation of the number of particles in the collision processes 
(notice that A n+ i is indirect characteristics of intermediate stages of encounters and 
collisions). This equality makes it possible to relate A™ +1 to A n+ i and thus A' n to 
A n+ \ . It is easy to see that unambiguously simplest form of the relationships is 



or, equivalently, 

A™ +1 {t, z, vi , v n \ v n+ x) = ~~t~3 [ A n+1 (t,z,v 1 ,...,v n ,v)dv . (3.2) 

J Ai{t,z,v)dv J 

This relationship does not touch on the correlations between members of the left-side 
n-particle cluster. 

Expression ( 13. 2\\ is a weakened version of the hypothesis of molecular chaos. It incor- 
porates the spatial statistical correlations of colliding particles, i.e. it asserts that only 
their velocities and momenta are statistically independent, but not their coordinates 
(thus, all the n + 1 particles may be mutually dependent in the configuration space). 
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Along with (13. 2p . Eqs J2.5l form a closed - we wish to stress this closure - hierarchy 
of kinetic equations. In the limit of spatially homogeneous gas, this hierarchy per- 



Moreover, this equation follows already from the first of of formulas (13.21) when we note 
that the conditions of inter-consistent normalization of the set of DF, 



(with Q being total, infinite, volume of the system), can be reduced in the homogeneous 
limit to the local form J F n+ idv n+ i = F n (F = 1). This implies J A n+ idv n+ i = A n , 
and from (13. 2p with n = 1 we have A™(t, v±\ v 2 ) = Ai(t, f 1 )A 1 (t, v 2 ) . 

In general non-homogeneous case, such relations no longer hold, since the exact 
global form of the conditions of mutual consistency of DF cannot be replaced by a 
spatially local form. The evolution of one-particle DF is determined by the entire 
infinite chain of equations (12.51) and (13.21) . and becomes definitely non-Marcovian, in 
contrast with evolution in the Boltzmann model. Clearly, this then leads to a low- 
frequency temporal dispersion of the spatially non-local kinetic transport coefficients 
of the gas. In turn, this dispersion may serve as a source of information about the 
low-frequency fluctuations of the kinetic coefficients, as we will see below. Of course, 
it would be wrong to think about the spatial non-homogeneity as a cause of these 
fluctuations (they take place also in homogenous and equilibrium states). In fact, the 
non-homogeneity gives only the means by which they manifest themselves in the one- 
time DF F n , due to dependencies of F n in non-homogeneous non-equilibrium states 
on the kinetic coefficients. Next, let us consider a simplest no n- homogeneous problem 
concerning self-diffusion of gas particles. 



mits the completely factored solution, A n (t, v 1 , v n ) = YYj=i Ai(£> v j) 5 an d becomes 
equivalent to the Boltzmann equation 
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4. 1/f noise of self-diffusion 



To analyze self-diffusion we need to eliminate from the kinetic equations the hydro- 
dynamic modes associated with the five integrals of motion of the system as whole. 
This can be done easily by taking the known formal approach (see e.g. [2]): replac- 
ing the probability distribution of the external particle in the collision integral by the 
equilibrium one-particle DF. In our notation, replacing (I3.2p by 



where A (v) = (27r)~ 3 / 2 exp (— 1> 2 /2) is the equilibrium Maxwell velocity distribution. 

Physically, this replacement describes a situation in which the gas is in equilibrium 
state in the macroscopic thermodynamical sense. There is only a small perturbation 
of the statistical equilibrium with regard to a single marked "test" particle and its 
immediate surroundings. The statistical state of the surroundings will be described by 
the set of DF which stand on the left sides of (I2.5P and (13.11) . The rest of the gas serves 
as the thermostat. Clearly, if the one-particle DF is assigned to the test particle, then 
the higher-order DF represent clusters consisting of the test particle and n — 1 other 
particles from its surroundings. Thus we assign the index 1 to our test particle. 

Further, we also will make use of the Green-Kubo theorem, according to which (see 
e.g. [6]) the generalized diffusion coefficient D(t, V) , which figures in the general 
non-local form of the self-diffusion equation (see e.g. [T2"]). 



- with V = d/dR and W(t, R) denoting probability density distribution of coordinate 
of the diffusing particle, - can be related to the linear response of DF of the test particle 
to an infinitely weak (potential) force f e x{<li) applied to it (we take in mind, as usually, 
that the force is "switched on" at some time moment, e.g. t = , before which the gas 
was in all respects at the equilibrium). Namely, it is a straightforward matter to prove 



A™ +1 {t,z,vi , ...,v n \v n+1 ) = A (t,v n+1 ) / A n+1 (t, z,vx,...,v n ,v)dv , 



(4.1) 




(4.2) 
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that the following relation holds (see Appendix A): 
where f ex {k) is the Fourier transform of f e x{qi) and 

/•oo 

D(p,V) = / rfr5(r,V) . 
./o 

To find the response we have to add terms m~ l f ex (q\) dF n / dv\ to the left sides 
of the BBGKY equations. In the kinetic equations (13.11) . i.e. in the framework of 
the collisional approximation under the low-density limit (and in the dimensionless 
notation), these terms look as f ex (z) dA n / dv± . The replacement of q\ by z in the 
argument of f ex presupposes that a change of f ex over length scales < a is negligibly 
small and that the external force has no substantial effect on the dynamics of collisions. 
These assumptions are well legitimate in the limits af ex /T — > and /U — > if, in 
addition, characteristic spatial scale (wave length) of the force, k , is not too small: 

i^r 1 > Ao. 

After the relations (14.11) are substituted into Eqs J3.ll the Boltzmann collision integral 
on the right sides transforms into a generalized Fokker-Planck operator A (or the 
"Boltzmann-Lorentz operator" [2]) defined by 

^j*P( v j) = J dv n+1 S jn+1 ip(vj) Ao(v n+1 ) . 
As the result, we come to the equations (in the dimensionless form) 

(4.4) 



d 



3=1 J V 3=1 J 



with the equilibrium initial conditions 

n 

A n \ t=0 = A° n = Y[A ( Vj ) . 

3=1 

In principle, D can be found directly from evolution of an non-equilibrium state in 
absence of external forces. This approach, however, is less convenient since it requires 
special consideration of initial stage of the evolution before its kinetic stage. 
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Next, consider the mentioned linear response, setting on A n = + n ( — > ), 
taking the Fourier transform in z and the Laplace transform in t , and denoting Fourier 
transforms of <p n by n . Then Eqs J4.4l yield 

P H I 0n = 5Z A J / ^ n+1 dtVl+1 + u i/exP _1 ^° • (4.5) 

n j=l ) 3=1 J 

Since these equations cannot be solved in their general form, we will simplify the 
problem. First, we restrict the analysis to the first two terms in the expansion of D 
in the gradient of the force inhomogeneity: 

D(p,ik) = D (p) + {ikfD x {p) + ... . 

Second, we choose the simplified model form of the operator A , namely, the Einstein- 
Fokker-Planck (EFP) operator: 

Here we have used that Ao(vj) = . 

In reality, such the operator could not arise from a short-range interaction potential 
(since it corresponds to scattering through only small angles). This circumstance, how- 
ever, should not corrupt the qualitative side of our results, because all that is required 
of A is that it ensures relaxation of velocities to thermal equilibrium. This choice is 
convenient in that all the eigenf unctions of A turn to product of Aq and polynomials. 
We may recall that the operator A corresponding to the Maxwell interaction potential 
has the same property [2j (therefore the following calculations can be generalized to 
this case). But the EFP operator has a further advantage: it makes it possible to 
exploit the separation of variables and thus to work with only projections of vector 
variables onto the force's wave vector k , i.e. to deal with formally one-dimensional 
problem (notice that if f ex is a potential force then vector f ex \\k). 
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Our next step is expansion of the response in the wave vector of the force-induced 
inhomogeneity: 

oo 
N=0 

Using formula (14. 3 p to relate this series to the above expansion of the diffusion coeffi- 
cient, we find 

A) = J viCf ] dv x , D x + Dl/p = J V1 C[ 2) d Vl = 5 1 . 
Inserting the same series into (14.51) . we come to equations 

n „ 

PC^ ] = / C%dv n+1 + Vl A° n , 

j=i J 

n „ 1 n 

PCS" = E A W C { nldv n+l - C^Wvi • 

i=i J i=i 
The equations from the first row can be solved easily with our choice of A and yield 

C™=^, A,(P)=-^-- (4-6) 
p + 7 p + 7 

The solution of the second row of equations for C$p should be sought in the form 



= \ + /3 n j A° n D n , 

where a n and (3 n are functions of p alone. For them we find the equations 
p(3 n = 27« n+ i , pa n = - 2'ja n+1 - l/n , 

and then 

1^1/ 2^y- n a i 



a. 



P ~~ n J \ V J Pn 



Now consider the functions 

5 n = / v dv\ ... dv r , 



the first of which determines D\ (see above). Multiplying the chain of equations for 

(2) 

C n by v\ and integrating over all velocities, we find, after some algebra, 
pS n = -7<W - n' 1 [ (n + 2) a n + p n ] D . 
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Substituting the expressions found above for a n and (3 n into this equation, we find 
8\ and thus D\ in the form of the repeating sum: 



D 2 
V 



p2 



n=l 



11- 



_L 1 00 1 



71 A — ' 7 

J =n J 



where X = 7/p. Then, with the help of the identity 



r x 

n - l X n ~ l = X- 1 / y n - 1 dy 
Jo 



we can transform the series over n to the easy summable form: 

n rX 00 



D 2 
V 



/ Ec-vr 1 

J n=l 



p 2 X 

Do r /p fin (1 + 2/) 
P7 7o 



?2 



1 



ln(l + y) 



+ 



1 + y 



(4.7) 



2/ l + 2y 

We have gone into the details of the calculations to demonstrate the characteristic 
p dependence of the response which is rather general and holds under another choice 
of the operator A . Let us consider this dependence at \p\ <C 7 and thus examine 
the behavior of the diffusion coefficient at low frequencies. From (14.61) and (14. 7p at 
p/7 — > we have D — I/7 and 



D(p,ik) « D 



l + (^) 2 ^(ln 2 ^ + c) + ... 



(4.8) 



2p V P 

where c is a numerical constant. In the dimensional notation, evidently, we must write 
Dq = ^o r o/7 = (T/m)r m and replace 7 by 7/V = l/r m in the logarithm, with r m 
being the momentum relaxation time. 

Further, we turn to direct statistical characteristics of random diffusive ("Brownian") 
displacement of the test particle as described by the diffusion equation (14.21) . Knowing 
the first N terms of the expansion of the diffusion coefficient in ik , we can find 
theoretically the first N + 2 of the statistical moments of the displacement, 



M n (t) 



R n W{t,R)dR . 
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From (14. 2 p we have, particularly, 

DM 
p _ 

Now, combining these formulas with (14.71) and (14.81) and performing the inverse Laplace 
transform for t r m , , we obtain 



M 2 (t)e- pt dt = , [°° M 4 (t)e- pt dt = — 

Jo P Jo P 



Di{p) 



\ n 2± + c >l n ± + c » 

Tm. Tm. 



(4.9) 



M 2 (t) « 2D t , M 4 (t) « 3M 2 (t) + 6D^ 2 

where D = Tr m / m , and c' and c" are numerical constants. 

Let us compare this result with that which would follow from the canonic approxi- 
mation based on the molecular chaos hypothesis and the Boltzmann-Lorentz equation. 
At that, in place of our infinite system of equations (14. 5 p we would have the single 
closed equation, 

(p + ikv^fa = Ai0i + v 1 f ex p~ 1 A (v 1 ) , 
which at arbitrary choice of the operator A implies 

M 2 (t) w 2D t , Af 4 (t) w 3M 2 (t) + c A^ t , (4.10) 

with Co being a numerical constant ( cq = when A is of the EFP type) and the same 
diffusivity D as above (by its sense, D = lim t -+aa M 2 (t) /2t ) . 

The second term in the expressions for M&(t) in ( 14.91) and ( 14.101) is the fourth-order 
cumulant (semi-invariant) of the test particle displacement, x 4 (t) = M 4 (t) — 3M|(t) . 
As is known, it is a measure of "non-Gaussianity" of the displacement. In particular, 
it shows how substantially a (more or less random) value of diffusivity D measured 
from a concrete realization of the Brownian motion can differ from the diffusivity Dq 
characterizing the average over statistical ensemble of realizations. The "Boltzmann" 
asymptotic ( I4.10p . that is x 4 (t) oc t , means that random wandering of the test particle 
can be exhaustively described by the single statistical parameter D . 

This is not the case, however, if x 4 (t) oc t v with v > 1 . Such kind of the asymptotic 
means that it is no longer possible to pack an arbitrary typical realization within the 
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framework of a single parameter D . It is not difficult to make sure that such the 
asymptotic is statistically equivalent to the existence if flicker fluctuations of diffusivity 
with low-frequency power spectrum oc cu - ^ -1 ) [TJ El El QUI Ell- 
in the theory under consideration, in contrast with the Boltzmann theory, we are 
just in such the situation, since, according to (14.91) . oc t 2 In 2 (t / T m ) . It can be 

described, on a coarse enough time scale (when relatively small terms proportional to 
t in x 4 (t) can be neglected), as a Gaussian random walk with a random diffusivity 
D = D(t) . Treating D as global characteristics of the entire time interval (0,t) 
under observation, we come to a "doubly random" walk whose path probability density 
distribution and statistical moments can be written as follow, 

W(t,R) = ( [47vtD(t)}- 1/2 exp [-R 2 /4tD(t)\ ) , 
M 2 (t) = (2tD(t)) , M 4 (t) = 3([2tD(t)} 2 } , 

where the angle brackets mean average in respect to D with (D) = D . The corre- 
sponding interpretation of the asymptotic (14. 9p in terms of power spectrum, Sd{&) , 
of the diffusivity fluctuations at u>r m 1 yields (see Appendix B) 

D 2 S d (uj) « - In— . (4.11) 

Thus, the random diffusivity, or the "rate of diffusion" , of the test particle (and hence 
of any particle) is a random process with a 1/f spectrum. Let us discuss this result^ . 

A measurement of spectrum of the diffusivity fluctuations is nothing but a measure- 
ment of the equilibrium average value of definite fourth-order polynomial functional 
of the particle velocity. An experiment of this type is not merely thinkable but really 



3 The approach suggested twenty years ago in the present paper [ Sov. Phys. JETP 67 (12), 2469 
(1988)] later was simplified and improved in [ Yu. E. Kuzovlev, "On statistics and 1/f noise of Brownian 
motion in Boltzmann-Grad gas and finite gas on torus. I. Infinite gas", arXiv: cond-mat/0609515]. 
There the asymptotic 

Wit f ex P [-i* 2 /4Dt] w (D_)d_D = £(5/2) 

was found, with w(x) — x~ 3 exp (— 1/x) representing effective probability density distribution of 
the relative rate of diffusion x — D/D ■ 
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has been carried out by Voss and Clarke (see references in However, it is vastly 

simpler to imagine measurements of usual quadratic functionals of the drift velocity 
of gas particle (under influence of an external force) and thus spectra of fluctuations 
in particle's mobility. Natural reasonings prompt that these fluctuations should be (at 
least at low frequencies) a statistical copy of the fluctuations in diffusivity. In the Ap- 
pendix C we prove that this is indeed the case for the system under our consideration 
(the "extension of the Einstein relation to fluctuations" takes place). 

From the standpoint of a principal experimental test of the theory, it would be inter- 
esting to extend it to multi-component gases, in particular, weakly ionized gas. It can 
be shown that diffusivity fluctuations of some component of a mixture represent a gen- 
eralized flicker noise with power spectrum uj~ a , where the exponent takes a value from 
interval 1 < a < 2 , depending on ratios of masses and momentum relaxation times of 
the components. Hence, the same principal mechanism of diffusivity fluctuations may 
produces a variety of flicker-type spectra. This can be illustrated also be spectrum 
uj~ 1 ln~ 2 (1/cjto) found in the earlier phenomenological theory (see IU El EH])- 

The spectrum (14. lip , as well as the asymptotic (14.91) . give clear evidences that the 
diffusivity and mobility fluctuations behave like statistically non-stationary and non- 
ergodic random processes. This statement means that in measurements of the diffu- 
sivity or mobility on a concrete particular realization of the particle's random walk the 
variance (in the sense of the ensemble average) of the result does not decrease, or even 
increases instead, when increasing length of the observation. The probabilistic aspects 
of such a behavior were studied in P, [TT] . Here, it is necessary to emphasize that such 
the statistical non-stationarity has no relation to a ther mo dynamical non-stationarity, 
since formulas (14. 8 p and (14.91) concern thermodynamically equilibrium gas. The matter 
is that the non-stationarity manifests itself only in dependencies on duration of obser- 
vation of a test particle, but not on time moment when the observation starts. Thus, 
the diffusivity and mobility, as well as kinetic and transport coefficients in general, 
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can undergo a non-stationary low-frequency fluctuations even in thermodynamically 
stationary and equilibrium systems. This principal possibility should be taken in mind 
when interpreting experiments. 

5. Conclusion 

The 1/f self-diffusion noise found above is the property of an infinitely low-density 
gas. Hence, the mechanism of this noise by its very nature is indifferent to the gas 
density and is not related to any dynamic many-particle correlations (e.g. through 
repeated collisions). It can be related only to fluctuations in the rate and efficiency 
of collisions of a given gas particle with the rest of gas because of randomness of 
geometrical factors of their encounters. The corresponding noise is just the 1/f noise 
if the system constantly forgets a total number of previous collisions of the particle 
and their past rates [9l dO, E] • More rigorously, the loss of memory is implied by the 
property of mixing of phase trajectories of the system in its full phase space. For a gas 
this property was proved already in [3]. Notice that in fact in [1] it was shown also 
that generally just the mixing of phase trajectories causes their non-ergodic behavior, so 
that time-average rate of collisions (or some other events) along a particular trajectory 
is not obligated to coincide with the ensemble-average rate, regardless of duration of 
the time averaging. The essential paradoxical point is that although such the random 
behavior results from the loss of memory, it can be described in the statistical language 
only by means of infinitely long-living correlations. 

In the derivation of kinetic equations (like the Boltzmann equation), however, always 
an assumption is made (like the molecular chaos) in order to replace actual random 
rate of collisions, or other elementary kinetic events, by some ensemble-average valued 



This real mission of the molecular chaos hypothesis (Boltzmann's "Stosszahlansatz" ) still is not 
realized even by important physicists. The common naive opinion is that it claims independence of 
particles' velocities. But in fact it introduces also a priori predetermined rate of collisions (although 
neither the underlying mechanics nor the statistical ensemble if initial conditions to this mechanics 
do present such a quantity). In contrast to it, the "weakened molecular chaos" ansatz, as expressed 
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[9l [10]. This decisive assumption is thus introduced already into the "zero-order" 
theory of infinitely low-density gas. Of course, taking into account finite gas density 
effects (see e.g. [3]) results in appearance of characteristic long-scale hydrodynamic 
fluctuations which contribute to both velocity and diffusivity of any particle and thus 
to inter-particle correlations. But, in spite if these complications, the above pointed 
out mechanism of the 1/f fluctuations anyway remains lostL. 

In the present work, our goal was just a correct exposure and reconstruction of the 
loss, at least in the low-density limit, at that accenting principal role of spatial non- 
homogeneity of gas. Notice that the generalization of the genuine Boltzmann equation 
to inhomogeneous situations by means of automatic adding of the drift term already has 
attracted critical comments more than once in the past (see e.g. [5]). The hierarchy 
of kinetic equations found in Sec. 3 is nothing but the correct (compatible with the 
concept of collision) formulation of the Boltzmann equation for inhomogeneous case. 

At last, let us note a relation of our theory to the formally exact generalized non- 
Marcovian kinetic equation for one-particle distribution function obtained in [6J by the 
method of projection operators. In [5] it was shown that analyticity of the kernel of 
this equation in respect to the Laplace transform variable (p in our notation) means 
its coincidence in the low-density limit with the Boltzmann equation. Hence the latter 
is invalid if the kernel is nonanalytic. The asymptotic (14. 5 p indicates just such the case. 
Thus, our specific theory agrees with the abstract theory of [6]. 



above in Sec. 3 by the equality A™ +1 (t, z,vi,...,v n \ v n+ i) = A\ (t, z, v n +i)A' n (t, z,v\ ,...,v n ) , does not 
dictate a certain rate of collisions, since generally A' n (t, z, V\ , v n ) ^ A n (t, z, V\ , v n ) . 

5 This fact prompts that, analogously, the kinetic theory of electron-phonon systems in solids (in 
particular, the theory of polarons) loses 1/f fluctuations in the phonon relaxation rates and electron 
mobilities at very beginning of the theory, when it attracts the hypothesis of "random phases" or 
another ansatz to throw out statistical inter-(quasi-)particle correlations. About that, see e.g. [8], 
the article [Yu. Kuzovlev, "Relaxation and 1/f noise in phonon systems", JETP 84 (6), 1138 (1997)] 
and preprint [Yu. E. Kuzovlev, "Kinetic theory beyond conventional approximations and 1/f noise", 
arXiv: |cond-m at / 9903350j] . Thus, in essence, the present paper (along with pj [11]) shows the way 
to rather general explanation of the 1/f- noise phenomenon. 
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Appendix A. 

The switching on at time t — a potential force field f ex ( r ) — —dU(r)/dr acting 
onto the test particle is described by the Hamiltonian 

H(t) = H o + U( qi )0(t) = H + 9(t)J2 l -^^-e ikqi = #o-]Tx fc (t)Q fc , 

k k 

where 9(t) is the Heaviside step function, H is Hamiltonian of the unperturbed 
gas, and Q k = e lkgi and X k (t) = — [ifc / ecc ( A;)] A; - 2 0~ 1 (9(t) in couples play roles of 
generalized variables and conjugated external forces. Let us consider the "flows" 

= ik Vl e ikqi , 



dQ k 



* ~ dt 

where v\ is velocity of the text particle. According to the Green-Kubo relations, to 
the first order in the forces we have 

(/*(*)> = f (h(t),h(t')) x k ,(t')dt' , 

where (...)o means the average over equilibrium statistical ensemble of phase trajec- 
tories of the system corresponding to the Gibbs canonic statistical ensemble of initial 
states of these trajectories, and (...) means the non-equilibrium average. In our case, 
in view of the homogeneity and isotropy of the equilibrium gas, the above relation takes 
the form 

v 1 e- ik ^F 1 (t,q 1 ,v 1 )^-dv 1 = (^e^ 1 ) = (A.l) 
1 ^ (v 1 (t)exp{- i k\q 1 (t)-q 1 (t')]}v 1 (t')) dt l 1 




The distribution of probability density of the test particle displacement, or path, in 
equilibrium gas can be determined by the expression 

W(t-t ,R) = (S( qi (t)- qi (t )-R)) . 
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Let us apply to it the differentiation operation d 2 /dtdto and then make the Fourier 
transform in respect to R (at to = 0) and the Laplace transform in respect to t . 
Using (vi )o = , we find 

POO 

p{pW(p,ik)-l] = -k / dte- pt (v 1 (t)ex P {-ik[q 1 (t)-q 1 (0)}}v 1 (0)) k , (A.2) 

Jo 

where W(p,ik) is the transform of W(t,R) . For it from (14. 2 p we have the identity 

pW(p, ik) — 1 = —k D(p,ik) kW(p,ik) , 

which, as combined with ( 1A.2I) . yields 

fn/^L , 2 = 1 I™ dte-^^t) ex V {-ik[ qi (t)- qi (0)]} Vl (0) ) . (A.3) 
p + D(p,ik) k 2 p J 

At last, performing in ( 1A.1I) the Laplace transform and comparing the result with ( 1A.3I) . 
we come to the formula (14. 3p . 

Appendix B. 

The asymptotic of the fourth-order cumulant *c±(t) expressed by (14.91) indicates non- 
stationary character of the fluctuations in diffusivity (observed on a particular random 
walk of the test particle). In other words, this asymptotic indicates the absence of a 
certain time-average value of the diffusivity. Therefore, D(t) should be treated as a 
non- stationary random process (in the sense of the probability theory). Since at t < r m 
we have M2 n (t) oc t 2n , this process begins from zero value: -D(O) = 0. The power 
spectrum of a random process of this sort is characterized, as is known, by so-called 
structural function: 

POO 

{ [D(t) - D(0)] 2 ) = 2 / [1 - costut] S d (lu) 

Jo 

At t ^> r m , after differentiation in respect to t , this expression and formulas of Sec. 4 
together yield: 

A 00 n / x • du d MAt) Dl , t 
/ ujS d (uj) smut — = — —— - w — In— . 
J 7T at 12 t z t T m 

The spectrum (14.111) follows immediately. 



duo 
2^ 



Appendix C. 
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Let at time t = a constant force / begins to act upon the test particle ("con- 
stant" means that / is independent on t and qi), and W(t,R;f) denotes resulting 
probability distribution of the particle's displacement, R(t) = qi(t) — ?i(0) . According 
to the generalized fluctuation-dissipation relations (see e.g. [13] and corresponding use 
of these relations in [S]), we can writer 

W{t, R- f) exp (-fR/T) = W(t, -R; f) . (C.l) 

It is not difficult to transform this exact equality into 

W(t, R; f) - W(t, -R; f) = tanh (§) [ W(t, R- f) + W(t, -R; f) ] . (C.2) 



Multiplying this by R and integrating over R, we find such a consequence of ( 1C.1I) : 

fR(t) \ 



R(t) ) = ( R(t) tanh ■ 



2T / 

Now, consider the third derivative of this equality with respect to / at / = : 



(C.3) 



\d 3 (R(t)y 


3 


~d 2 (R 2 (t)y 






[ dp 


f=o ZJ 


[ dP 




4T 3 



(C.4) 



and apply this general relation to our problem. The average displacement 



( R(t) ) 



dt' 



under the influence of the constant force f ex = f = const is determined by the homo- 
geneous solution to equations ( 14.4D . In the homogeneous case, however, they evidently 
reduce to the ordinary Boltzmann-Lorentz equation. We can thus assert that the aver- 
age drift velocity, ( Vi(t) ) , reaches saturation at t > r m . Consequently, the left side of 
( 1C4j) increases linearly with the time. However, the last term in (IC.4I) . which contains 
the fourth-order statistical moment of equilibrium displacement, (R A (t) )o = M^t) , 



6 Besides [13], the generalized fluctuation-dissipation relations were investigated e.g. in the works 
[G. N. Bochkov and Yu. E. Kuzovlev, "Fluctuation-dissipation relations for nonequilibrium processes 
in open systems", Sov. Phys. JETP 49, 543 (1979); "Fluctuation-dissipation theory of nonlinear 
viscosity", Sov. Phys. JETP 52. No. 12 (1980); "Nonlinear fluctuation-dissipation relations and sto- 
chastic models in nonequilibrium thermodynamics. I. Generalized fluctuation-dissipation theorem" , 
Physica A 106, 443-480 (1981)]. 
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due to (14. 9p increases far more rapidly. Hence, this term is compensated (exactly at 
t/r m — > oo ) by the first right-hand term of (1C.4I) . As the result, from (1C.3I) and (1C.4I) 
we obtain 

(R(t)) = ^M 2 (t) « ^/t , 

( R\t) ) = M 2 (t) + M 4 (t) » 2D t + ( 5 2 (t) ) . (C.5) 

On the other hand, we can write 

(R(t)) = t{fz(t))f, (R 2 (t)) = M 2 (t) + t 2 (^(t))f 2 , 

where //(£) is the mobility referred to the entire observation interval as a whole. A 
comparison with flC.5j) shows that the Einstein relation between the diffusivity and the 
mobility holds not only in the ordinary sense (for the ensemble-average values) but also 
for their fluctuations. Correspondingly, the spectrum (14.111) simultaneously refers to 
the relative fluctuations in mobility. A similar relation, concerning the spectral density 
of equilibrium electric voltage noise and electrical conductivity, has been confirmed 
experimentally in the famous Voss-Clarke experiment (see e.g. [8]). 
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